Planar approximation for spin transfer systems 
with application to tilted polarizer devices 
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Planar spin-transfer devices with dominating easy-plane anisotropy can be described by an ef- 
fective one-dimensional equation for the in-plane angle. Such a description provides an intuitive 
qualitative understanding of the magnetic dynamics. We give a detailed derivation of the effective 
planar equation and use it to describe magnetic switching in devices with tilted polarizer. 
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I. INTRODUCTION 

The spin-transfer effect is a non-equilibrium interac- 
tion that arises when a current of electrons flows through 
a non-coUinear magnetic texture^^3 ^ Spin transfer torque 
can lead to current induced magnetic switching in multi- 
layer devices or domain wall motion in devices with con- 
tinuous change of magnetization. Both phenomena serve 
as an underlying mechanism for a number of suggested 
memory and logic applications. 

Magnetic dynamics in spin-transfer devices can be de- 
scribed by the Landau-Lifshitz-Gilbert (LLG) equation. 
Analytic solutions of LLG can be easily found in the sim- 
plest case of easy axis magnetic anisotropy. However, 
when the form of anisotropy energy becomes more com- 
plicated the investigations of the stability of static equi- 
libria become much more involved. A study of the pre- 
cession cycles is even more complicated and often makes 
it necessary to resort to numeric simulations. Due to the 
complexity of the LLG equation it is always interesting 
to consider cases where some simplifying approximations 
can be made. 

In many devices the easy plane anisotropy energy is 
much larger than the other anisotropy energies, and the 
system is in the planar spintronic device regime^ (Fig.[l]). 
The limit of dominating easy plane energy is character- 
ized by a simplification of the dynamic equations^ , which 
comes not from the high symmetry of the problem, but 
from the existence of a small parameter: the ratio of the 
energy modulation within the plane to the easy plane en- 
ergy. The deviation of the magnetization from the plane 
becomes small, making the motion effectively one dimen- 
sional. As a result, an effective description in terms of 
just one azimuthal angle becomes possible. 

Our publications^^~ introduced the planar approxima- 
tion in the presence of spin-transfer torques and listed 
a number of results highlighting the practical use of the 
method. In this paper we give a detailed derivation the 
effective planar equation for a macrospin free layer in the 
presence of spin transfer torques (Sec. lIIII) . We then show 
in Sec. IIVI how this equation can be applied to a system 
with the tilted polarizer and obtain a qualitative picture 
of the device dynamics. 



II. MAGNETIC DYNAMICS OF THE FREE 
LAYER 

We consider a conventional spin-transfer device con- 
sisting of a a magnetic polarizer (fixed layer) and a small 
magnet (free layer) with electric current flowing from one 
to another (Fig. [T]) . The free layer is influenced by the 
spin transfer torque, while the polarizer is too large to feel 
it. It is assumed that due to the large exchange stiffness 
the free layers can be described by a macrospin model. Its 
magnetic dynamics is described by the Landau-Lifshitz- 
Gilbert equation with the spin transfer torque term. 

The state of the free layer is characterized by just one 
vector, its total magnetic moment M = Mn which has 
a constant absolute value M and a direction given by a 
unit vector n{t). The LLG equation^ii^ reads: 



Sn 



- u{n)[n X [s X n]] + a[n x h] . (1) 



Here the re-scaled energy e = {■^/M)E has the dimen- 
sions of frequency and is expressed through the total 
magnetic energy E{x\) of the free layer; 7 is the gyro- 
magnetic ratio, and a is the Gilbert damping constant. 
The second term on the right is the spin transfer torque. 
Unit vector s points along the direction of the polarizer 
moment, and the spin transfer strength u(n) is propor- 
tional to the electric current /— . In general, spin transfer 
strength is a function of the angle between the polarizer 
and the free layer u(n) — g[{n ■ s)] /. The spin current 
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FIG. 1: Planar spin-transfer devices. Hashed parts of the 
devices are ferromagnetic, clear parts are made from a non- 
magnetic metal. 
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efRciency factor ^[(n • s)] is a material and device specific 
function. 

In polar angles {9, (f) equation ^ reads 
1 de 

9 + a(t)sm0 - 



)SinW — ab 



sin 6 
de 
09 



+ u(n)(s • ee) , 
+ u{n){s-e^), (2) 



where the unit tangent vectors eg and are defined in 
the Appendix 1X1 

We assume that the easy plane is defined by 9 = 7r/2, 
and the rescaled magnetic energy has the form 

e = ^Jcos2 6' + £r(6',0) ■ 

The first term represents the easy plane anisotropy, with 
the planar frequency ujp related to the easy plane con- 
stant Kp as ujp = jKp/M. The remainder is the 
"residual" energy. The planar limit is characterized by 
Wp — oo. Large easy plane constant forces the en- 
ergy minima to be very close to the easy plane and 
the low energy solutions of LLG to have the property 
e{t) = 7r/2 -I- S9 with 69 -> 0. Equations © can then be 
expanded in small parameters 

M«i and 

OJp UJ„ 



Up 



By truncating this expansion one obtains the effective 
planar approximation. 



III. DERIVATION OF THE EFFECTIVE 
PLANAR EQUATION 



Next, we make the approximations. In the Wp — ^ oo 
limit the solution is expected to have a property 9{t) = 
TT/2 + d9{t) with 69 0. Expanding all quantities on the 
r.h.s. of (0 in small 66 up to the first order we get 



69 



fg" + Fg^6e - a{ujp5e + f ° + Fl69) 
I + a2 

tjp69 + F^ + Fl69 -f a{F^ + F^69) 



1 + 

where we have used the notation 

7^0 



(6) 
(7) 



In the approximation (I6I7P equations are linear with re- 
spect to the unknown function 9{t), but still fully non- 
linear with respect to 4'{t). 

Eq. ([7]) can be viewed as an equation on 59 which gives 



69 = 



(I + a^ 



Fl - aFS 



ujp + Fl + aFl 



(8) 



so that the out-of-plane deviation becomes a "slave" of 
the in-plane motion.'^ The presence of the large ujp in the 
denominator ensures the smallness of 69. Substituting 
the resulting expression 69 = q{4>, <P) back into equation 
([6]) one obtains a second order differential equation for a 
single unknown function (pit) 



dq ■■ dq 



F" ~ aF" 



aojr, 



aFl 



Explicitly separating the large easy plane terms, we 
rewrite equation © as 

9 + acj) sin 9 = /e + ug 

(j)sm9 — a9 = — cjpCos^sin^ + + 

where the residual anisotropy is responsible for the terms 

I dsr 



smb 
'89 ' 



(3) 



and the spin transfer torque produces the terms 

ue{9,4>) = u(n)(s • eg) , 

u^{9,4>) = u(n)(s • e^) . (4) 

We also introduce a notation Fg^^f, — /e,^ + ug^^, and 
re-write the LLG system as 



Fg — a(— cjp cos 9 sin 9 + F^) 

— Wp cos 9 sm9 -\- Frj, + aFg 
(1 -I- a^) sin 6* 



(5) 



By denoting il{(j)) — uip + F^ + aFg and simplifying the 
terms we get 



1 + 



+ a 



= F^- 



dq 



(1 + a' 



n 



F^iF°,+aF^) 



(9) 



In this form the equation is still too complicated to be 
useful but since it was obtained from the approximations 
(|6l7p we are allowed to drop the terms that are smaller 
or equal to the already neglected ones. Those were the 
F"69^ terms in the expansion of F, and the third order 
terms ujp69^ in the easy plane energy expansion. While 
formally the latter are of higher order in the 69 expansion, 
the large coefficient LUp can causes them to have the same 
order of magnitude. 

To compare the orders of magnitude of the terms con- 
sistently, we need to know the order of magnitude of 69. 
At the present stage we know that 66 is small but its 
exact order of magnitude is not known because we do 
not have an estimate for the (f> term in the numerator of 
Eq. ©. 



A. Simple residual energy in the absence of spin 
torque 



B. Arbitrary residual energy in the absence of spin 
torque 



To estimate 59, let us first consider the problem in the 
absence of spin transfer {F — f),^ assuming a simple 
form of residual energy = £r{4>)- In this case we find 
fl = fl^n^O and get 



(l + a2)0-a/o 



UJv 



(10) 



Equation ^ takes a form 



ujp \ ujp dtp J 

Assuming that the residual energy £^(0) does not have 
any special points of fast change we can estimate Fg ^ 
e 'ti LOp. Then 

ujp ocp 

and we can approximate the equation by 



1 



OJr, 



a<i) = le = — TTT 



The equation above has the form of the Newton's equa- 
tion for a particle of mass (1 + Q!^)/wp moving in a one- 
dimensional potential Er (</>), subject to a viscous friction 
force with a friction coefficient a. Our goal is to estimate 
the value of 0, i.e., the speed of the "effective particle". 

The particle's characteristic speed depends on the total 
energy and on the relative strength of the friction forces. 
We will assume that the energy is of the order of Sr (this 
is the mathematical equivalent of our original assumption 
about the low-energy dynamics of the moment). Further- 
more, in the present paper we will concentrate on the case 
of a — >■ that corresponds to an almost frictionless mo- 
tion of the particle. Then one can use the approximate 
energy conservation and write {l/2ijjp)(fP' = Sr for the 
maximum speed. This gives 

Using similar arguments one can estimate the maximum 
acceleration as 

Note that the viscous friction can be approximately ne- 
glected when a0 <C Fq ^ Sr, i.e., the Gilbert damping 
constant a has to be not just small compared to unity 
but satisfy a more stringent inequality 



Let us return to the approximation ()6I7|) with the gen- 
eral form of the residual energy Sr = er{6,4>) and zero 
current, u — 0, F — f. It is now possible to use the 
a posteriori estimate (|12p for 69 to consider the orders 
of magnitude of the terms on the right hand sides of 
the equations. We start the discussion from the "slave" 
equation ([7]). Here 



af « erJ — 



af^S9 < Er 




Up 



As we see, the orders of magnitude of the terms form a 
series 



er_ 

UJ,-, 




(13) 



with the general term given by Srisr/i^p)^^'^- 

The terms neglected in transition from system ([5]) to 
d?]) were 

ujJ9^ - 



56' ^ Er — 



This means that in ([7]) one should only keep the terms 
of the order Er and higher. Lower order terms would 
be comparable to some of the discarded terms. Using 
this argument we discard f^59, af^ and af^56. The 
1/(1 -t- Q!^) factor in ^ can be expanded using Eq. ([TT|) 

Er 



= 1 -I- ,5, 5 - < — 

1 + a'' Up 

This inequality shows that 1/ (l + a^) can can be approxi- 
mated by unity in Eq. ([7|) without changing the accuracy. 
After those simplifications equation ([8]) takes the form 



a < 



< 1 . 



(11) 



We see now that </> is the largest term in the denominator 
of and as a result obtain an estimate 



(12) 




As for the equation (jG)), the terms discarded in going 
from dS]) to ([6]) were 



aujpd9^ <^ Er — 

UJp 
Er 



5V 
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IS now 



and therefore we have to keep the terms of the orders 
£r\/ ErjiOp and higher. Thus the and a/° terms 

should be kept in (jS]) but the aF^66 terms should be 
discarded. One can also conclude that it is safe to re- 
place the 1/(1 -l-a^) factor by unity. Equation 
replaced by 

where one should use 

j-n 



UJp 



since the term fgf^/ij^p ^ s^/ujp is of the same order as 
the already discarded terms. 

Re-deriving Eq. ^ in this approximation one gets 



1 



'4. IS 



9e^(7r/2,(^) 



This is the equation for the effective particle as discussed 
in the previous section, except that the viscous friction 
coefficient seems to acquire a correction. While the cor- 
rection is of the order Sr/wp, it is added to a small num- 
ber a and thus can potentially have a significant effect of 
changing the sign of the dissipation term. However, one 
finds 



d dsr d 



de,. 







(14) 



so the correction actually vanishes. We come back to the 
effective equation 



acf) ■ 



9e,-(7r/2, i 
del) 



which corresponds to the most natural generalization of 
the equation derived in the previous section for a spe- 
cial form of the residual energy Sr — Sr{(f>)- The posi- 
tive effective friction coefficient ensures that the effective 
particle always stops in the point of energy minimum, as 
expected for a closed system with dissipation described 
by the LLG equation without spin torques. 



C. Effective equation in the presence of spin torque 

Finally, we proceed to the derivation of the effective 
equation in the presence of the spin torque. Consider 
approximations (|6l7p with Sr = £r{9, 4>) and u ^ 0. 

The order of magnitude of the extra terms produced 
by the current will depend on the value of u. One con- 
dition that certainly has to be satisfied is the smallness 
of the spin torque compared to the anisotropy torques 
produced by the easy plane contribution. The latter are 
responsible for the terms of the order y^evw^ in Eqs (j6l7l) . 
Thus it seems that u should not exceed Sr, which is the 



largest term before y/ErUJp in the series (|13p . Such a 
conclusion is correct for a general situation. We will, 
however, see below that in some special cases the current 
can be increased up to it ~ ^er(^p without violating the 
dominance of the easy plane anisotropy torque. 

To include those cases we assume u ^ ^Sr'^p and 
revisit Eqs. (I6I7I) discarding the terms smaller than 
er\J Crl'^-p in Eq. ([B]), and smaller than in Eq. ([7]). 

Eq. ([7]) which acquires the form 

^ = (c^p + u\)8d + ]l + u%. 

(as in the previous section, one can prove that the factor 
\|(\^r o?') can be approximated by unity without loss of 
accuracy). By solving for 8B and expanding the denom- 
inator up to the same accuracy we find the form of the 
slave condition ([8]) 




(15) 



Differentiating both sides one gets 







Up J 


LOp 








' del) 



1 d[ulul]\ 



del) 



- (16) 



Returning to Eq. (jG]) we find that with the declared ac- 
curacy it can be rewritten as 

Substituting 59 from (jlSp and discarding any terms that 
are smaller than £r^j£r/^p, we get 



fl+ul , ulul 



11^1 + ulJl + ulul , ululul 



p 



The last step is to use Eq. (|16p to express 69 on the left 
hand side. This gives the form of the effective equation 
([9|) without the terms below our accuracy 



"0 \ 



p j 




ulu\ 



Hi! 



Je + I 5~ 
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where identity ([T4| can be used to simplify the bracketed 
expression on the second Une. 

We now cast the effective planar equation in its final 
form 

• ^'^^ = -(^e//)' . (17) 



m0 + ae//(; 



Here primes denote differentiation with respect to 4>, and 
the parameters are given by 




^eff 



a 



(18) 



'TT 



U 



— U — Ug \ . 

UJp 

Equations PT)) and ((T5| constitute the first main result 
of this paper. 

In the presence of the current, u ^ 0, the right hand 
side of pT)) contains additional "effective forces" added 
to the —Er' term. Since all functions depend on just one 
variable 0, these forces can be always represented as the 
derivatives of an additional energy U according to the 
definition (fTS]). 

One of the forces, namely the Ug term, requires a spe- 
cial discussion. When u ^Je^ZJ^ this term becomes 
larger than Sr on the right hand side of Eq. ((TT)) . and 
the estimates for and made in Sec. IIII Al become in- 
valid. As it was discussed above, this means that in a 
general case with non-zero the effective equations ([TTl 
[T5)) can be only used for currents u< Sr- However, if Ug 
is identically equal to zero, while the other spin torque 
terms in U and aeff remain non-zero, one can apply 
Eqs. ([T71 [T5)) for currents up to u ~ ^ErUip- 

Corrections to the friction coefficient are explicitly de- 
pendent on the current magnitude. In the presence 
of spin torque the sign of the friction coefficient may 
change^"— , reflecting the possible influx of the energy 
from the current source into the system. 

Below we investigate the application of the effective 
planar equation to the device with a "tilted polarizer" 
geometry. 



IV. TILTED POLARIZER DEVICE 

In tilted polarizer devices vector s is pointed at an 
angle dg to the axis z and constitutes an angle i: /2 — Og 
with the easy plane. We will assume that s is in the {x, z) 
plane, i.e., 0s = 0. 

To calculate ug and from Eq. ^ one needs to know 
the function u{n). In many cases^'-^^ it has a form 

"(n) = -; 7 r 

l + 5i(n-s) 



We will consider the case of small gi ^ 1 and approxi- 
mate 

zi(n) =5o/(l-5i(n-s)) (19) 

Using the expressions in Appendix \^ we flnd 

ue — gol [I — gi{sin 6s sin 9 cos (j) + cos 6s cos 9)] 
X (sin 9s cos 6 cos cj) — cos 9s sin 6) , 
= — go^ [1 — 3i(sin6's sin^cost/) -I- cos6's cos^)] 
X sin 6s sin . 

Therefore 

Ug — —gal [1 — gi sin 6s cos 0] cos 9s , 

Mi = —gol [1 — gisin6s cos(f)]sin6sSin(f>. (20) 



and 



ul — — gQl[sm6sC0S(t> + gi{sm^ 9sC0S^ cos^ 6s)] 



-gol 91 sin 6s cos 9s sin < 



A. In-plane polarizer 



(21) 



In the case of in-plane polarizer with 9s — ix jl further 
simplifications happen: 



u% = 

,0 



u"^ = -(1 - 01 cose/)) sin 

v}g = — (I — COS 0) COS 



As we see, the in-plane polarizer is one of the special 
cases with Ug — 0, discussed at the end of Sec. IIII CI 
Consequently, the effective equation can be used up to 
the currents u ^ gol ~ ^ErtOp. The coefficients (fT8)) 
acquire the form 

I 

TO = , 

UJp 

.go/[2cos(/)-.gi(3cos2 0- 1)] 
"e// = a H , (22) 



-U' 



_ go/(l - gi coscf)) ^ 

UJp 

X [go/(l — 5i COS0) sin(/>cos(/) - 
-(/e sin 0-1-/^ cos ( 



(23) 



Importantly, the (p^ term in Eq. (1171) vanishes identically. 

In Refs. 6,7 the in-plane polarizer was considered in 
the case of gi = and residual energy 



sm f cos 



h sin 9 cos i 



(24) 



describing a device with small easy axis anisotropy uja <C 
UJp in an external magnetic field /i, both pointed along the 
X axis. In this case one finds /g = / ° = and expressions 
reproduce the results obtained in Refs.0 and[3- 
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B. General case of a tilted polarizer 

When the polarizer magnetization s points at an arbi- 
trary angle 6^ the term is nonzero and we have to limit 
the current magnitudes to g^)! ^ Sr to maintain the valid- 
ity of Eq. (fT7|) . With smaller currents more terms can be 
discarded from the effective equation without changing 
its accuracy. Parameter expressions reduce to 



m 
-U' 



1 



K)' + ui 



(25) 



Moreover, for qqI < the nonlinear term with (j)^ be- 
comes small enough to be dropped from Eq. pT| . 
The effective planar equation now reads 



de. 



eft 



(26) 



where the effective damping and the effective energy can 
be expressed through the polarizer tilting angle using 
Eqs. (I2ni), (EH), and ^ 

aef f = a + (2 sin 6s cos (j) ~ gi[3 sin^ 9s cos^ </> — 1]) , 



Eeff = er{-,(t>) + goIicOsBs 



gi sm Us cos t)s sm ( 



(27) 



C. Switching diagram of the tilted polarizer device 

Let us now discuss the consequences of the modifica- 
tion of a — >■ aef f and Sr — > £e/ / in the presence of spin 
torque. The advantage of the effective planar approxi- 
mation is the possibility of using the analogy with the 
particle motion which enables one to use mechanical in- 
tuition to qualitatively predict the behavior of the solu- 
tions of Eq. (|26p , and thus understand the dynamics of 
the spin-transfer device. 

We will assume the standard nanopillar device de- 
scribed by the residual energy (|24p . In the special case 
of an in-plane polarizer this problem was discussed in 
our earlier publications j^iii^ Modifications of the effec- 
tive damping transform the particle motion qualitatively 
when Ueffifj)) changes sign. Modifications oisef /{((>) be- 
come qualitatively important when equilibrium points 
appear or disappear as the energy profile is deformed. 

Our goal here is to generalize the results of Rcfs. J..J,,9 
to the case of a non-zero polarizer tilt and show that 
the effective planar approach allows one to understand 
the qualitative picture of the motion without doing the 
detailed calculations. The case of a tilted polarizer was 
considered in recent publication o^^'^'^ using conventional 




FIG. 2: Effective energy profile in the case of a tilted polarizer. 



methods. The switching diagrams were calculated as- 
suming that s is tilted straight up from the easy axis 
direction. This corresponds to our assumption of 0s = 0. 
It was also assumed that there is no magnetic field, h = 0. 
We adopt the assumptions of Refs. [H and[l3| to illustrate 
the application of the formalism developed in this paper. 

In the case of a general tilt angle 9s the terms con- 
taining a small factor gi give negligible corrections in the 
expressions (|27p for tte// and £eff- They can only be- 
come important when 9s approaches zero or tt/2 and the 
main terms vanish. We will assume that s is not to close 
to either the in-plane or the perpendicular directions and 
an inequality cos 0^, sin 0s 3> gi holds. Then we can use 
the simplified form of Eqs. (|27p 

2 go / sin ^s 
aef f = a-\ — ■ cos cp , 



Eeff = er{-,(f>) + golcos9s ■ (f) 



(28) 



Finding the critical currents in the narrow bands of an- 
gles 0s ~ or w 7r/2 where the gi terms are important 
would require a more careful consideration. 

The profile of the energy eeff{(f>) is shown in Fig. [2] 
The spin torque produces the tilted washboard potential 
in which the effective particle moves. The washboard 
tilt reflects the fact that spin torques can change the 
magnetic energy of the system by using work performed 
by the current source. The energy minima are shifted 
from their zero current positions at = (parallel, or P 
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Ml 




P 



FIG. 3: Effective particle performing a 360° motion, starting 
from the energy maximum Mi and reacfiing an equivalent 
maximum M3. The red (gray) interval denotes the region of 
Oe// < 0. 



(a) Ml (b) Ml 




AP AP 



FIG. 4: (a) Effective particle performing small oscillations 
(IPP precession) near the minimum corresponding to the AP 
state that lies within he red (gray) interval of Ofe// < 0. 
(b) With growing current the amplitude of the oscillations 
increases and the particle reaches the M2 maximum point. 
Above this threshold the IPP precession is unstable. 



state) and 4> = tt, (antiparallel, or AP state) to the new 
positions (j)o{I) that are given by the equation 



sin 200 



2,90^ cos ( 



One can check that the energy minima located at the an- 
gles 4'Qmm + 7™ are still separated by the energy maxima 
located at at ipomax = 7r/2 - 0omjn + 7rn (Fig. [2b)- To 
be concise will still call the (jjomin minimum the P-point, 
and the (j)omin + ti" minimum the AP point. 

As the current grows, the washboard tilts more and 
more, until the extrema or eeff{(t>) disappear altogether 
(Fig. [2t). A short calculations shows that this occurs at 
a critical current 



2go cos 9s 



(29) 



For / > /»£ the effective particle slides down the slope 
of the potential energy profile regardless of the sign and 
magnitude of cteff- This corresponds to n performing 
full 360° rotations in the azimuthal angle 0. In the spin 
transfer literature such a regime is called the out-of-plane 
(OPP) precession. 

Importantly, the OPP precession can exist even at 
/</»£. When the particle moves down the washboard, 
the drop of its potential energy during one spatial pe- 
riod may be large enough to overcome the friction energy 
loss (Fig. . Therefore there must be a second critical 
current Iqpp < I*e, such that the OPP precession hap- 
pens for / > Iqpp- In the interval Iqpp < I < I*e the 
stationary equilibrium of the particle at the energy min- 
imum coexists with the state of OPP precession. The 
functional form of Iopp{Ss) depends on the energy pro- 
file and the friction coefficient. Our goal here is not to 
find the expression for it, but to see how far can we pro- 
ceed in qualitative understanding of the device dynamics 
without doing the actual calculations. 

The effective damping tte// changes its sign in the 
vicinity of the angles 04.„ = 27r(n -I- 1/2) for J > and 



4>-n = ^im for / < when the current magnitude ex- 
ceeds the threshold 



2(7o sin 6 s 



(30) 



The regions of Ofe// < are shown in Figs. I2|3I bv red 
(gray) thicker line. On the intervals of negative friction 
the system acquires energy instead of loosing it. This 
makes it easier for the particle to achieve the state of 
the OPP precession. The actual calculation of the Iqpp 
threshold must take into account the energy gain due to 
the tilt of the potential and the presence of negative fric- 
tion intervals. Both features mathematically represent 
the ability of the spin torques to transfer energy from 
the current source to the system. 

Another important effect of negative aeff is the lo- 
cal destabilization of the energy minima. If the mini- 
mum point 4>a{I) lies within the interval of Oieff < 0, 
it becomes unstable and small oscillations around it are 
developed. Fig. 2^ shows such oscillations near the AP 
minimum which, according to Eq. (j28p . is destabilized for 
/ > 0. These oscillations correspond to the precession of 
n around the equilibrium point and are called an in-plane 
(IPP) precession in the spin transfer literature. 

The amplitude of the oscillations is determined by the 
balance of energy influx and dissipation on the intervals 
of negative and positive friction/^ii'^ As the current is 
increased, the amplitude grows and eventually becomes 
so large that the particle reaches the crest of the poten- 
tial (Fig. |3]d) and falls down into the neighboring valley. 
This process leads to the destruction of the IPP state. 
The latter therefore exists between the two threshold cur- 
rents. For the AP equilibrium these are the Iap thresh- 
old, where the AP point becomes unstable, and the Ijpp 
threshold, where the precession amplitude becomes too 
large to be contained in the AP valley. 

The critical current of minimum point destabilization 
is determined from the equation agf f{(t>oiI)) = which 
can be rewritten as 



Iap cos[(I)o (Iap)] 

/pCOs[0o(/_p)] = 



(31) 
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FIG. 5: Critical lines of the P/AP equilibrium states and the 
IPP precession state on the {1, 6s) switching diagrams. Each 
line is marked by the name of the state which gets destabilized 
on it. The name is put on the side of the line where the state 
is stable. Inset: an impossible arrangement of the AP and 
IPP lines. 



for P and AP minima. The equations show that the 
threshold currents satisfy |/| > /*q in both cases. This 
result can be naturally understood as follows. The fric- 
tion first becomes negative at the (j) — or cj) — t: points 
at / = ±/*Q,. But the P and AP minima are shifted from 
the 0, TT points to the (j>o{I) points. In order to destabi- 
lize them the negative friction interval has to grow large 
enough to cover the actual minima positions. 

The critical current // pp depends on the shape of the 
potential in the entire interval traveled by the particle in 
Fig. Its actual calculation is not the goal of our qual- 
itative approach. We can make two general statements 
about 1 1 pp. First, the destabilization of the precession 
state certainly happens at a current that is larger than 
the one required for the destabilization of the correspond- 
ing energy minimum. For example, Ijpp > Iap- 

To set the stage for the second observation, we proceed 
to the discussion of the switching diagram. In our case 
the experimental parameters are the current / and the 
tilting angle 0s. Various critical currents are represented 
as lines on the {I,ds) plane and divide it into domains 
with different sets of stable states. The currents 
Iap, and Iapp are sketched as functions of the angle 
9s in Fig. [5] The lines I^ei^s) and /*q(0s), given by 
Eqs. ([^ and ([50]) . intersect at a certain point. Due to 
Iap > I*a the Iap{Os) has to be located to the right of 
the I^a{Ss) line on the diagram. It intersects the I^:^{9s) 
line as well — one can prove that Iap{Os) is a decreasing 
function of the tilt angle. The Iipp{Os) > Iap{Ss) is also 
a decreasing function and the corresponding line has to 
be located even further to the right. Both Iap{Os) and 
Iipp{Os) cross the I*e{ds) line, and one can prove that 
they do it at the same point T. This is second general 
property of the Ijpp threshold. 

The second property can be proven by considering a 
hypothetical switching diagram shown in the inset in 
Fig. [5] where it is assumed that that the AP and IPP 
lines cross at different points Ti and T2. Consider 

the point T2. Approaching it from within the domain of 



existence of the IPP precession one should observe a 
decreasing amplitude of oscillations around the AP min- 
imum, because the size of the valley around AP shrinks 
to zero. At T2 the amplitude of sustained oscillations 
around AP should be equal to zero. At the same time 
the AP point should be unstable because the current is 
larger than the Iap threshold. But an instability of a 
minimum with zero amplitude of oscillations is the def- 
inition of the Iap current, so T2 has to lie on the AP 
line. This argument proves that the assumption about 
the existence of two different crossing points Ti^2 was in- 
consistent. 

The relationships between Iap and Iipp discussed 
above follow from the fact that both currents are de- 
termined by the energy and friction near the same min- 
imum. The lopp current depends on the details of £e// 
and aef / on the whole 27r interval of the angle (j) and thus 
no qualitative relationships for it can be found, except 
for the already mentioned Iqpp < I*e{Os)- A qualitative 
sketch of the full switching diagram is given in Fig. [51 
The / < part of the sketch can be obtained by reflect- 
ing it with respect to the vertical axis: the diagram is 
symmetric with respect to the / — > — / transformation. 
This is a consequence of two symmetries built into the 
energy and friction functions (|28l) : the 7r-periodicity of 
the energy ([M)) . and the ae//(7r — (j),—I) = aeff{4>,I) 
symmetry of the effective friction. The latter depends on 
s being directed in the symmetry plane and the fact that 
the gi terms were dropped. If either of the two symme- 
tries is violated, the P and AP states would no longer be 
equivalent in all respects. 

Notably, Fig. [5] reproduces all qualitative features of 
the switching diagrams obtained in Refs. [l^ and [l^ by 
conventional methods. In addition, it brings important 
qualitative understanding of the behavior of critical cur- 




FIG. 6: Switching diagram of a spin-transfer device with a 
tilted polarizer: a qualitative sketch. Each line is marked 
by the name of the state which gets destabilized on it. The 
arrows next to the name point to the side of the line where 
this state is stable. Gray horizonthal strips denote the regions 
where the approximation (|28|l may fail and (|27p should be 
used. 
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rents as a function of system parameters and approxi- 
mations used for the spin torque efficiency factor 5(11). 
For example, we see that the threshold currents in the 
two most frequently considered limiting cases of the per- 
pendicular (0s — 0) and in-plane (Og — tt/2) polarizers 
will be most sensitive to the ^(n) function used for the 
calculations. 

Finally, we should mention that the results of the ef- 
fective planar approach are not limited to the qualita- 
tive conclusions discussed in this paper. It allows one to 
calculate the critical currents, often being the only one 
providing analytical expressions in the case of precession 
states. In our case the /*e threshold is given by the ex- 
pression ([29|). The Ip and Iap currents can be obtained 
from Eqs. ([3l|). At the crossing point 9s* of the AP, IPP 
and /*£ lines both equations should hold, which gives 



tan — 



In the limit of small friction a ^ ^Jujajujp considered 
here the critical currents Iopp{Ss) and Iipp{9s) related 
to the precession states may be found using the method 
introduced in Ref. l2|. 



merging and disappearance of the minimum and max- 
imum points of £eff- For Og close to 7r/2 the destabi- 
lization happens locally, changing the nature of the AP 
equilibrium from stable to unstable. This type of com- 
petition is not unique to the systems with strong easy 
plane anisotropy — it was shown in Ref. [3 that it may 
happen in any spin-transfer device. 
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FIG. 7: Definitions of the tangent vectors and polar angles. 



V. CONCLUSIONS 

We have given the detailed derivation of the effective 
planar equation for spin-transfer devices with dominat- 
ing easy plane anisotropy and illustrated its application 
by performing a qualitative study of a spin-transfer de- 
vice with tilted polarizer. Once the parameters of the 
effective equation are found, the approach allows one to 
understand the dynamics qualitatively without perform- 
ing detailed calculations. This is especially important in 
the case of precession cycles which are usually studied 
numerically. The method also elucidates the role of ap- 
proximations used to model the spin torque and shows 
the limits of their applicability. 

The obtained switching diagram demonstrates a com- 
petition between the two types of switching. For small 9s 
the destabilization of the AP minimum results from the 



Appendix A: Vector definitions 

We use the standard definitions of polar coordinates 
and tangent vectors (see Fig.[71): 

n = (sin 9 cos 0, sin 9 sin 0, cos 9) 
ee = (cos cos 0, cos sin — sin 6') (Al) 
= (— sin0, cos (/), 0) 

For the polarizer unit vector s with polar angles {9s-, 0s) 
the scalar product expressions are 



(s • eg) = sin6's cos0cos( 
(s • e^) = sm9sSm{(f)s - 



cost/, smf 
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